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Abstract 

The master equation of quantum optical density operator is transformed to the equation of characteristic 
function. The parametric amplification and amplitude damping as well as the phase damping are considered. 
The solution for the most general initial quantum state is obtained for parametric amplification and amplitude 
f~*) ' damping. The purity of one mode Gaussian system and the entanglement of two mode Gaussian system are 

fv^ , studied. 
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(^ ; 1 Introduction 



■^^ ' Quantum information with continuous variables (CV) [T] [5] is a flourishing field, as shown by the spectacular 
^ . implementations of deterministic teleportation schemes [3] [1] [S] [H], quantum key distribution protocols ^Tj, 
entanglement swapping [S] [5], dense coding [H], quantum state storage ^U\ and quantum computation [TT] 
processes in quantum optical settings . The crucial resource enabling a better-than-classical manipulation and 
^+ ' processing of information is CV entanglement. In all such practical instances the information and entangle- 
f^ , ment contained in a given quantum state of the system, so precious for the realization of any specific task, is 
^>0 ' constantly threatened by the unavoidable interaction with the environment. Such an interaction entangles the 
system of interest with the environment, causing some amount of information to be scattered and lost in the 
f^ ' environment. The overall process, corresponding to a non unitary evolution of the system, is commonly referred 
Mh. to as decoherence [4] [5] . In this work we study the decoherences of general states of continuous variable systems 
whose evolutions are ruled by optical master equations. We will consider the parametric amplification, ampli- 
tude damping and phase damping [12j of a general quantum CV state in the fashion of quantum characteristic 
^ ' function. The main starting point for this research work was the result of Lindblad's bounded generator of 
^^, a completely positive quantum dynamical semigroup |13j . Quantum characteristic function had been used to 
treat the amplitude damping of one mode squeezed states [14j . 
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vh ■ 2 Time evolution of characteristic function 

The density matrix obeys the following master equation [12] [13] [E] 



with the quadratic Hamiltonian 



| = 4[i^,p] + (A+/:.)p. (1) 

H = ?i^ -(%fea]4 - V*jk"-jak) (2) 

jk 

where 77 is a complex symmetric matrix. In the single-mode case, this Hamiltonian describes two-photon 
downconversion from an undepleted (classical) punip[T5]. The full multi-mode model describes quasi-particle 
excitation in a BEC within the Bogoliubov approximation T6]. This item represents the parametric amplifier. 
While the amplitude damping is described by £1, 



Ap = ^^{(n, + l)L[a,]p + n,L[a]]p (3) 

3 

-w*M[aj\p ~ WjM[a]]p}. 

where the Lindblad super-operators are defined as -/j[o|p = 25/90' — o^op — po'd and il/[o|p = 2dpd — oop — poo. 
The requirement of positivity of the density operator imposes the constraint \wj\ < nj(nj + 1). At thermal 
equilibrium, i.e. Wj — 0, rij is equal to the average thermal photon number of the environment. If Wj ^ 0, then 
the bath j is said to be 'squeezed' fl7|. 
£2 describes the phase damping, 

^2P = E Y^[«l«.^]/^- (4) 

j 

We now transform the density operator master equation to the diffusion equation of the characteristic func- 
tion. We use characteristic function because it is more convinent for our problem. One may use Glauber's 
P-representation instead, but Glauber's P-representation does not exist for some of the states (e.g.[14] [E]), 
although the generalised positive P-representation does always exist [19'. Any quantum state can be equiv- 
alently specified by its characteristic function. Every operator A S B{T-C) is completely determined by its 
characteristic function XA '■= trlADip)] [20 , where I?(/^) = cxp(/ia^ — p*a) is the displacement operator, with 
H = [//i, /j,2, • • • , Ps] ,a— [oi, 02, • • • , fls]-^ and the total number of modes is s. It follows that A may be written in 
terms of XA as [21]: A = /[Hi ~~^]XA{i^)'D{—p). The density matrix p can be expressed with its characteristic 
function x- X = ^''[p^(m)] ■ Multiplying V{p) to the master equation then taking trace, the master equation of 
density operator will be transformed to the diffusion eqation of the characteristic function. It should be noted 
that the complex parameters pj are not a function of time, thus -^ — ir[-g|I?(/i)], the parametric amplification 
part in the form of characteristic function will be [TS] 

dx , * dx 



jk jk ^^ ^^^ 

The master equation can be transformed to the diffusion equation of the characteristic function, it is 



dx sr^f * 9x , * dx . ,„^ 

Tt - -\^'W^^^^+ri,,H^) (6) 
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2z^-u-'a|p, 



p,| -^ + ((271, + 1) |^,f - w*p) - w,pf)x} 
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Where we denote pj as \pj\e^^\ and we should take care about that the variables are pj and p* in the 
amplification while they are \pj\ ,6j in the damping. 

3 Solutions of some special cases 

Firstly, let us consider Tj — jj — for all j, the solution to the amplification is 

X{p,p*,t) =x{'^,i'*,0), (7) 

with 

/ cosh* (hi i) _^*£i£M^M \ 

V ^r^V cosh(|77|i) J 



where the matrix cosh and sinh functions are defined as|22j 

"iih|«l = ' + ^«' + l{«'f + ---. (9) 

Then suppose 77 = and "fj = for all j, the solution to the amplitude damping equation of x is 

X{^i,t) = x(A^e-S^,0)exp{-^(l - e-i^^*)[(n, + i) \fi/ - ^w*fi^ - ^w.fif]}. (10) 

j 

Where /J,e 2 is the abbreviation of (/iie 2 , ^26 2 , • • • , /ZgC 2 ). Next, suppose 77 = and Fj = for all j, 
the solution to the phase equation of x then will be 



^(A.,Ai*,i)=y'dxx(A^e",/i*e-",0)n(27r7.i)"'/'exp(-^), 



(11) 



where /xe" is the abbreviation of (/iie*^i,/i2e*^^, • • • ,/iie*^=). The simultaneous amplitude and phase damping 
(77 — 0) for any initial characteristic function is (for Wj = 0) [23] 

x(a^,M*,0 - /'dxx(Me-*+",M*e-*-",0)n(2^7,0"'/' (12) 

exph^-(l-e-i^'*)(n, + i)|M,n. 
The density matrix then can be obtained by making use of operator integral. 

4 The parametric amplifier and the ampUtude damping 

The diffusion equation ^ now is with its "fj = for all j. Suppose the solution to the diffusion equation is 

x(M,M*,i)=x(^,^*,0)exp[i(z.,-z.*)(^^ J ) (7.*, -z.f - i(^, -M*) (^ ^ j)(^*,-^f], (13) 

where ly — fiM + ij,*N, with M and N being time varying matrices, a and f3 are constant matrices and a^ — a, 
(3 — 0^ . M and N are the solutions of the following matrix equations 

^ - -,fN-\M. ,», 

^ = -IM-^N. (15) 
where T = diag{Ti^ r2, • • • , Fg}. While a and [3 are the solution of the following matrix equations 

2770 + 2a*r; - F/3 - /3F + Fw + uiF = 0, (16) 

Fa + aF-27;*/3-2/3*77-F(7i+i)- (7T+i)F = 0. (17) 

where w = diag{'Wi ,W2, ■ ■ ■ ,Ws},n — diag{ni , 772 , • • • ,77s}- The constant matrices a and /? can be worked out 
as the solution of linear algebraic equations p^ and ^T7\ . What left is to solve matrix equations p^ and p3|) . 
There are two situations that the equations are solvable. The first case is that all the modes undergo the same 



amplitude damping, that is Fi = r2 
and (fTSll have solution 



M ^ 

N = 



Fs, thus F = Fil^. F commutes with any matrix. Equations p4p 
-i^Uosh*{Mt), (18) 



ir,sinh(|77|^ 



1^1 



(19) 



The solution to the one-mode situation is simple. In the two-mode situation, M and A'' can be further simplified. 
For 7? is a symmetric matrix, we can express it as rj = r/ouo + ''71 ci -|- rj^as, where ctq = I2, cr^ (i = 1, 2, 3) are 
Pauli matrices. The a2 term is nulled by the symmetry of t]. By using of the algebra of Pauli matrices, we arrive 
at the results: 



cosli(|?7| t) 

smhilrjl t) 
— n — '/ 



(To cosh{^/{C + A)/2t) cosh{y/{C-A)/2t) 

+ sinh(v/(C + v4)/2t) sinh(v/(C-v4)/2t)"CT • ~b 
1 



-ao[sinh(V(C + B)t)/^y{C + B) + sinh( ^(C - B)t)/^{C - B)]r] 



(20) 
(21) 



1, 



-|--[sinh(V(C + B)t)/^{C + B)- sinh(V(C - S)i)/ ^(C - B)]l? ■ & 7^ 



Where C = jr^ol -I- |?/i| +1%! , A = h'/o ~ '7i ~ ^^l L -S = \/C'^ ~ A"^ ■ fo is a unit vector and it is equal to 
(^^i + '7i'7oi *%'?! ~ ^ViVh VqVs + VsVo)/ ^- Thus the characteristic function is simplified. If the initial state is 
Gaussian. The correlation matrix of the time dependent state can be obtained explicitly. 

The second case is that 7y is a real matrix while the amplitude damping can be different for each mode. The 
solution of equations (HH) and p^)) will be 



M = -[exp(-?7t-y)+exp(?7t-y)], 
1 Tt Tt 

N = -[cXp(-77t-y)-CXp(7;t-y)]. 



For two-mode situation, this M and N can be simplified to 

1 



M 



N 



-e '^^*[cosh{Bit)ao - smh{Bit)'a ■ bi] 



1 
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cosh(_B20'''o ~ sinh(i320 cr • &2], 



r 
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[cosh(i?if)cro ^ sinh(i3ii) a ■ hi] 



(22) 
(23) 

(24) 
(25) 



--e"'^"*[cosh(B2t)(To - sinh(B2i)"CT • 62 



where Ci,2 = ±m + Wi + r2); ^1,2 = ^rjl + {\{Ti-T2) ±VzY] "?i,2 = (±r?i, 0, ±773 + i(ri - F2))/Si, 



5 One mode Gaussian system 

In this and next sections, we will apply the solutions of previous section to Gaussian states. In the case of the 
parametric amplifier and the amplitude damping, we consider the case of Fi = F2 = • • • = Fg = F (hereafter F 
is a number, not the matrix used in previous sections). If the initial state is Gaussian, its characteristic function 
has the form of x(a'jM*:0) = exp[/z77i'''(0) — /i*m-^(0) — ^(/i, — /i*)7(0)(/x*, — /z)-^], the state will keep to be a 
Gaussian state in later evolution. The complex correlation matrix (CM) 7(0) should be chosen in a fashion that 
the intial state is physical. The time evolution of the complex CM is 



lit) 



-,-rt 



cosh*(|77|i) ri — ^ 

sin\i(\ri\t) 
^1 ' 



+ sinh(|77|i) 

cosh(|77| t) 



(7(0) 



a 13* 
P a* 



^, sinh(|77|t) 

SiHiliJlM^ cosh(f^li) 



cosh*(|ry|t) rj ^ 



-t- 



a (3* 
(3 a* 



(26) 



The time evolution of the complex first moment m is 



.-irir ,^^„„„u*n„u^ , ™*/n^ si^Mhl i) ^ 



i(t)=e-^^*[m(0)cosh*(|??|t)+m*(0) ;" ' ij]. (27) 

\m 



The solution of the characteristic function for one mode system is characterized by Eq. (fT8|) and Eq. (fT9|l 
where rj is simply a complex number now, and (for F 7^ 2 jr;!) 

a - a* = _J_^[r(rT+i) + r;*z^ + ,7Z^*], (28) 

^ = ^, \, M Mn+h + ir^~2\rj\^)w + 2rj'w*]. (29) 

1 ^ — 4 I77I ^ 

The result reduces to that of [H] when w = 0,n = and ?7 is real. The degree of the mixedness of a quantum 
state p is characterized by means of the so called purity fip = Trp^ . With characteristic function, we have 

IJ,p — J ^-^ \x{tJ;P*)\ ■ For one mode Gaussian state, it reads 

" 71 (t) 72 (i) 



Pp{t) = , with j{t) 

2^(Re7i(i))^-|72(i)l' 



72 (t) 7i(i) 



(30) 



We now consider the evolution of coherent state, squeezed state and thermal state. The coherent state is given 

by exp[m(0)a^ — m*(0)a] |0) ,the complex CM is 7(0) = diag{^, ^}. The squeezed state is given by exp[(Ca^^ — 

(*a'^)/2] |0) , where ( = r exp[j(/?],the complex CM is 7(0) = ^[cosh(2r)tTo + sinh(2r) cosipai + sinh(2r) sm(pa2]- 

The thermal state is given by X]^o[^/(-^ ^ ^)]" 1*^) ('^1 '^ith N is the average photon number, the complex 

CM is 7(0) =diag{N+^,N+^}. 

For F > 2 r/l , all these states (and other Gaussian initial states) will tend to a Gaussian state with complex 

\ a B* ] 
CM 7(00) = „ after a sufficient large evolution time. The ultimate purity is 



When the phase angle of the amplification 77 is equal to the phase angle of the 'squeezed' environment w , the 
ultimate purity is maximized, which is 

Mpmax(oo) - \{ ^ A n+\f - |^|'}-i (32) 

^ 1 ^ — 4 |?7| ^ 

The ultimate purity is a monotonically decreasing function of the amplification \r]\. 

For F < 2 jryl , we consider the case of w = and t] is real and positive for simplicity. When the initial state 
is a thermal state with 7(0) = diag{N + ^, N + ^} (coherent state is the special case of iV = when CM is 
concerned). Denote coshro = 2r\l ^ ^jf — F^, then 

7(t) = e'^\N + -)[cosh(2r7i)o-o + sinh(2r;i)cri] 

+ (n + -)sinhro[e~'"*(sinh(277i + ro)o-o + cosh(277i + ro)cri) — (sinhroCTo + coshroCTi)]. (33) 
The state is a squeezed thermal state, with its purity 



where N' = N +\,v! = (n+ \) sinhro. For large i, ^p(t) « \ exp[-(77- F/2)i]/v^7^^+WiVV-^. The purity is 
an decreasing function of the amplification r\ at large t. When the initial state is a squeezed state with squeezed 
parameter C = rexp[i(/3], we consider the simple case of (/? = 0,then 

^p(i) ^ i{[ie(2')-nt+2r +^'gro(g(2^-r)t _ 1)] [ie-(2'?+n*-2r ^ ^/^-ro (^ „ ^_(2^+r)t )]|- i ^ (35) 



For large t, ppit) w \ exp[-(77 - F/2)t]/V7i'2 +7i'e2'~-''o/2. 



6 Two mode Gaussian system 



The algebra equation of a and (3 in two mode system is complicated in general situation. To investigate 
the entanglement property of the amplifier, we will set rjo — r]^ — which corresponds to no single mode 
amplification. Thus r) = j/iCTi, it corresponds to two mode amplification. The solution is simplified to (with 
Fi = r2 = r, and the two modes undergo the same noise (n — 7iol2) for simplicity) 



r2 



P 



4 hi I 
1 



r(no + i), rjiWb + iilwl 
rjiWa+rilwl r(no + ^) 



r(r2- 4^110 



(r2-2|,7inu;, 
2r;ir(no + i) 



m'>^b 



2r/ir(rTo + i 



Wb 



vlK 



(36) 
(37) 



where we have denoted w = diag{wa, Wb} for the two modes a and b. If the initial state is Gaussian, the state 
will keep to be a Gaussian state in later evolution. 

For F > 2 1 771 1, the state will tend to a Gaussian state which is characterized by the residue complex CM 
7(00) after a sufficient larger evolution time. The Peres-Horodecki criterion [24] [25] for separability of the state 
is an inequality on the real parameter CM 7^6 (cxd), which is 



7re(oo) 



7a 7c 
T 

7c 76 



L 



with 



L = ^ 



1 

3 



i 

f 

i 

1 



^ 

-f 





a P* 
(3 a* 



LK 



-1 



(38) 



(39) 



Thus we have ja = aa<Jo-Im(3aO-i+Ref3a(^3, jb = ab(Jo-Im(3bcri+Re(3b(J3, 7c 
where we have denoted 



aa 

< 


Oic 
Oib 


, /? = 


Ha 


Pc 
Pb 



=ReacO-o— Im/JcCTi —ihnacO'2 +Re/3c0'3 , 

(40) 

The separable criterion takes the form of det7adet7h + (i -\Aei"fc\f-tr{'-^aJlcJlbJl]^ J) > |(det7a + Aei"fb) 
[24], with J = i(72. Hence in the form of a and /3, it will be 

,1 



det7;det7^ + (- 



|det 7^1)2 - ir(7>37c<^376'^37c^CT3) > ^(det7; + det70. 



where 



ai P* 
Pi a* 



a.b, c. 



(41) 



(42) 



We now consider the situation of w = for simplicity. Then aa 



ab 



r^(no + i 



r2-4|r,i| 



T,ac 



0;/3c 






Pa = Pb — 0, the CM 7re(oo) Can be transformed to the standard form [25] by local rotation. The standard 
form CM is 7,^^(00) = p2 41 ^ri (rcro ^S ctq + 2 \r]i\ ai (Xjcts). The inseparability criterion reads aa — \Pc\ < 5, which 
is 

FrTo<hi|<F/2. (43) 

The possible entanglement appears only when no < ^ . The entanglement of formation (EoF) of the inseparable 
state will be [2S] [?7] 

Ef = g{A{2aa - 2 |/3,|)) = g{A{^^^^)), (44) 

where A(z) = (z + z~^ — 2)/4 and g{x) ^ [x + \) log2(a; + 1) — a;log2 x is the bosonic entropy function. g{x) 
is a monotonically increasing function of x, thus the entanglement is a monotonically increasing function of 
the amplification |77i|( note that aa ~ \Pc\ < \)- The entanglement tends to its supremum when no = and 
|r?i| -> F/2. The supremum is Ef^ = g{\) = 0.5662 . 




Figure 1: Entanglement of formation of the two mode amplification and damping, with the environment noise 
no = 0-4 • The initial state is vacuum state. 



For r < 2 |ryi|, we consider case of w = and rji is real and positive for simplicity. From equations (|20|1 and 
1]), the time evolution solution of complex CM will be 



7(0 



,-rt 



cosh(77it)(To <8i ctq + smh{riit)ai <g) ai]['y{0) - (aaCo ® (Tq + PcCi ® cri)] 



[cosh(77ii)cro o-Q + sinh(77it)(Ti (g) cti] + aa^o (Xi ctq + PcCFi ® cti. 



(45) 



Where Pc is real and positive. If the initial state is a two mode squeezed thermal state with real CM 7re(0) 
{N + i)(cosh2rcro (8) cto + sinh2rcri ® 0-3), then 7,-e(i) — 7rei(i)o'o ® o'o+ 7re2(i)o'i ® era, with 



(46) 
(47) 

(48) 
(49) 

If the initial state is vacuum state, then A^ = 0, r = 0. Inequality ((15|) reduces to aa — /3c < 5 -It can be further 
simplified to 



7rei(i) = e ^'[(iV+-)cosh2(7yit + r)-aaCOsh(2?7it)-/3cSinh(277ii)]+aa 
7re2(i) = e~"[(7V+-)sinh2(77ii + r)-aaSinh(277ii)-/3cCosh(277ii)]+/3e 
The inseparable criterion for the state is 7rei(i) — 7re2(i) < ^, that is 

(Af + i) exph(277i + r)< - 2r] + (aa - /3c){l - exp[-(27;i + r)<]} < ^. 
The EoF of the inseparable state will be 

Ef(t) = g(A(27„i(i) - 27„2(i))). 
When ^ -^ 00, EoF saturates at 



ii;/(oo)=g(A(- 



-[2no + l])). 



?7i > maxjnoF, -}. 



(51) 



EoFs for no = 0.4, 0.6 are shown in Fig.l and Fig. 2. We can see that EoF is a monotonically increasing function 
of the amplification 771, Ef{oo) has the same expression for 771 > -^ and 771 < 2' (see (|44p '). Ef{oo) is shown in 
Fig. 3 (It is for any Gaussian initial state, not only for vacuum initial state). The inseparable criterion inequalities 
(H5)) and I^Tj) can be combine 

r/i/F > no. (52) 




Figure 2: Entanglement of formation of the two mode amplification and damping, with the environment noise 
TTo = 0.6 . The initial state is vacuum state. 




Ti/r 



Figure 3: The ultimate entanglement of formation with respect to amplification damping ratio rji/T and envi- 
ronment noise no 



7 Conclusion 

The master equation of quantum continuous variable system is converted to the equation of quantum character- 
istic function. It turn out to be a hnear partial differential equation about the characteristic function. The time 
evolution solution can be obtained exactly for any initial quantum optical state in several cases. The solvable 
cases include (1) parametric amplification [ 22] , (2) amplitude damping with thermal or squeezed noise |17| , 
(3) simultaneously amplitude and phase damping together with thermal noise P5]. (4) simultaneous multi-mode 
parametric amplification and amplitude damping with thermal or squeezed noise when each mode undergoes the 
same strength of damping, (5) simultaneous multi-mode real parametric amplification and amplitude damping 
with thermal or squeezed noise. 

The applications to one mode and two mode Gaussian initial conditions are investigated. In one mode 
Gaussian system, the purity of the evolution state is monotonically decreasing with the amplification. In the 
situation of less amplification {\r]\ < r/2), the ultimate purity is maximized when the phase of the amplification 
r] matches with the phase of the 'squeezed' environment w . In two mode Gaussian system, the entanglement 
of formation monotonically increases with the two mode amplification. In the situation of less amplification 
(jr^il < r/2), the supremum of the entanglement of formation is given. In the situation of over amplification 
{\rii\ > r/2), after a sufficient large time, the entanglement of formation saturates. The ultimate entanglement 
of formation is given as a function of the amplification damping ratio and the noise. It is independent of the 
initial Gaussian state. The ultimate state is separable when the amplification damping ratio is greater than the 
noise. 

In real experiment, the parametric amplification and the damping may occur successively. Thus the con- 
catenate of above solutions will represent most of the optical evolution system. 
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